Anomalous latent heat in non-equilibrium phase transitions 
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We study first-order phase transitions in a two-temperature system, where due to the time-scale 
separation all the basic thermodynamical quantities (free energy, entropy, etc) are well-defined. The 
sign of the latent heat is found to be counterintuitive: it is positive when going from the phase 
where the temperatures and the entropy are higher to the one where these quantities are lower. The 
effect exists only out of equilibrium and requires conflicting interactions. It is displayed on a lattice 
gas model of ferromagnetically interacting spin- 1/2 particles. 

PACS numbers: 64.60.Cn, 05.70.Ln, 75.10Hk 



The theory of equilibrium phase transitions is an es- 
tablished field with known achievements in describing 
transformations of various states of the matter 0, U 
Concepts borrowed from this theory apply for some non- 
equilibrium transitions which get equilibrium features 
on the macroscopic scale, e.g., equilibrium universality 
classes for non-equilibrium transitions , transitions to 
the glassy state analyzed with help of effective temper- 
ature 0,13, etc. Less is known, however, about phase 
transition scenarios which are impossible in equilibrium. 

Here we study such a truly non-equilibrium phase tran- 
sition scenario realized in the steady state of a two- 
temperature system. This system admits a natural gen- 
eralization of the equilibrium statistical thermodynamics, 
i.e., the quantities like entropy, internal energy, free en- 
ergy are well defined, because the constituents of the sys- 
tem are in local equilibrium. In spite of that, the system 
shows a counterintuitive type of first-order phase transi- 
tion, where the latent heat is positive {anomalous latent 
heat) when transforming a high temperatures (higher en- 
tropy) phase to the low temperature (lower entropy) one, 
i.e., in the first phase the energy is larger. This is a non- 
equilibrium effect. It is well-known, even from the every- 
day physics, that for equilibrium transitions the latent 
heat is negative li]. This plays a crucial role in the heat 
balance of the Earth and in the weather formation, since 
nearly 70% of the energy transferred from the Earth's 
surface is due to the latent heat consumed during the va- 
porization at the surface and released during the vapor 
condensation in the atmosphere 

First we shall recall the thermodynamics of two- 
temperature systems with different time-scales IE 13 ■ 
In contrast to the usual equilibrium case, the anoma- 
lous latent heat is not forbidden here. We then work out 
a simple model of mean-field Ising ferromagnet demon- 
strating the sought effect. A necessary condition for its 
existence is the presence of conflicting interactions. 

Consider a pair of coupled stochastic variables s and 
/ with Hamiltonian H{s,f), which interact with differ- 
ent thermal baths at temperatures Ts and Tf = T, re- 
spectively (s and / can denote a set of variables). For 
T = Ts — the stationary probability distribution 

P(s, /) of the system is Gibbsian: P(s, /) oc e~'^^'-^'^\ 
For T Ts we can derive the stationary P(s, /) if the 



variables have different characteristic times: s is slow, 
while / is fast (adiabatic limit). This derivation together 
with corrections coming from a large, but finite time-scale 
difference was given in based on stochastic equations 
of motion. Here we recall the heuristics [3, IE IS US ■ On 
the times relevant for /, s is fixed, and the conditional 
probability P{f\s) is Gibbsian: 
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Z{s) = T17 
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where Tr/ is the sum over all values of /. The steady- 
state P{s) is found by noting that on the times rele- 
vant for s, T is already in the conditional steady state. 
Thus the force dsH{s, /) acting on s can be averaged 
over P(/|s): Tiy [d,H{sJ)Pif\s)] = dsF{s), where 
F{s) = —T\nZ{s) is the conditional free energy. The 
steady P(s) is Gibbsian with the Hamiltonian F{s): 
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Tr,e^''=-^(''\ (2) 



and the common probability is P(s, /) — P{s)P{f\s). 

This two-temperature situation admits a (generalized) 
thermodynamical description, because in the adiabatic 
limit both variables are in local equilibrium 0, H 0- 
(It has certain analogies with non-equilibrium thermody- 
namics of the glassy state proposed in 0.) The average 
energy of the system isU = Ttsj [P{s, f)H{s, /)], or 
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(3) 



where the derivative is taken for fixed n and we used 
Eqs.im21). 2 is a generating function, since an arbitrary 
average is found via the proper source term analogously 
to ||2J). The entropies of s and / are, respectively, 

Ss = -Tt,P{s)\iiP{s), (4) 
S = -Tr,P(s) [Tr/P(/|s)lnP(s|/)] . (5) 

Eq.Q) is the usual definition of entropy. Eq.© is the 
conditional entropy; it appears due to the adiabatic 
limit. The total entropy is the sum of partial ones 
5tot = -Tr,j[P(s,/)lnP(s,/)] =5 + 5,. Now the 
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steady distributions 121 are obtained when minimiz- 
ing the free energy of the slow motion 



-TslnZ. 



(6) 



To this end, note from (^EJ that the free energy and the 
entropies are expressed as 

:f = u-ts-t,s,, (7) 

S, = -dT,T\T, S^-dTT\T^, (8) 

generalizing the usual thermodynamical relations. 

Let us now study the first-order phase transitions in 
this approach. Consider two phases I and h. Assume 
that when Tg decreases below some critical value ri'^' (T 
is fixed), the phase I dominates, since its free energy is 
smaller: Ti < Th- For Tg > 7s the dominating phase 
is h: Ti > Th- At high (low) temperatures Tg, I [h) is 

(c) . 

metastable. Ti — Th at Tg = Tg , while other quantities 
change by jump. Denoting 6X = A; — Xh we get 

6U - T,5Ss + T5S = -T,dT^ [5T]\t ~ Tdr [5T] |t, , (9) 

where 5U is the latent heat of the transition. Note that 
in general both 5Ss and 5S are functions of Tg and T . 
In the vicinity of Ts'^\ ST = Ti — Th is an increasing 
function of Tg (for fixed T); thus 5Ss < 0. The sign 
of 5S is, however, left open and we cannot conclude that 
5U < 0. It may be even positive (anomalous latent heat), 
provided 5S > Q, that is provided the system moves from 
I to h when decreasing T for constant Tg. 

Things get different if we decrease both Tg and T, with 
constant n = T/T^. Since ST{Ts,nTs) — Ti ~ Th has to 
be an increasing function of Tg, we get -^ST = (Sj; + 
ndT)5T > 0, and then shows that SU < confirming 
the equilibrium result for Tg = T l]. 

Thus for Ts the anomalous latent heat is not for- 
bidden, provided the transition is not driven by a pro- 
portional change of both temperatures. Here is a model 
demonstrating this effect. Particles are located at the 
nodes of a lattice embedded into a particle reservoir with 
chemical potential /i = —a < 0. Each particle car- 
ries an Ising spin. At short distances the particles re- 
pel each other so that not more than one particle can 
occupy a singe node. The Ising spins interact ferromag- 
netically. This interaction is active between two nodes 
only if they both are occupied; s = {si = ±1}^^ and 
f = {fi = 0, are, respectively, the spins and the 

occupations of N nodes. We assume all nodes can inter- 
act with each other (mean-field) . The Hamiltonian reads 



(10) 



where J/ {2N) > is the ferromagnetic coupling constant 
with J = C(l), as required for the extensivity. 

We make two assumptions. 1) s and / couple, re- 
spectively, to a spin bath and lattice bath at different 
temperatures, Tg = J- and T = i. 2) s (/) is slow 



(fast), i.e., the relaxation time Tf of / is much shorter 
than the relaxation time Tg of s. These times are driven 
by the interactions with thermal baths, since H contains 
only commuting terms; note that in the quantum setting 
Si = Si^z and fi — ^(1 -l-<Ti,z), where Si^z and ct^^z are the 
third Pauli matrices. 

The above assumptions are motivated by NMR/ESR 
physics, where for nuclear or electronic spins the spin 
bath is realized by relatively weak dipole interactions and 
does indeed lead to a temperature different from the lat- 
tice one This spin temperature can be tuned by 
external fields or by spin cooling and plays an important 
role in analyzing experiments The relaxation time 
Ts on which the spin temperature is established is known 
as T2 time in NMR/ESR, and usually varies between 
lO^'^s and Is depending on the material. The spin and 
the lattice temperatures tend to equalize on the Ti time, 
which for many magnetic materials amounts to minutes 
or hours ,11.] . Estimating the relaxation time Tf of / via 
the rotation-vibration mechanism as 10^^^ — lO^^s 11], 
we see that our assumptions are based on t/ ^ ^ Ti . 

The model is a suitable candidate for displaying the 
anomalous latent heat, since the two terms in the RHS of 
(|10|l are in conflict: for two lined up spins, SiSk > 0, the 
term —jJ2<ik>fifk^i^i^ tends to increase occupations, 
fifk > 0, while aJ2ifi makes this increase costly. The 
equilibrium (Tg = T) version of the model is related to 
Blume-Capel model 0- Similar models were employed 
recently for describing solid mixtures 0] , inverse freezing 
(reentrance) phenomena 12] and glassy physics [T^ . 

Employing eE. ^ c/dme'^'^'^^'+'^-^'^S.^'/-, 

where c = we express from Eas.(|21[Tn| the parti- 

tion function Z of the model as 



Z = Tr< 



Tr. 



Dme 



where Vm = YYl=i J dma, and we assumed that n = 
^ is an integer (replica method); later on we shall make 
continuation to real n. We also defined (j){m) = 1 + 
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.fc=±i 



.(11) 



Z is calculated by the saddle-point method, where one 
searches the deepest minimum of as a function of 
•nia- This minimum is reached for rua = m (replica 
symmetry), where m satisfies the stationarity condition 

r)m •^|m„=m = 0, i.e.. 



0"-i(m)e''"'™ - 0"-i(-m)e-''"''" 



0"(m) + 0"(-m) 
The corresponding value of T obtained from (|1H) is 



T 



Tgln [0"(TO)-H(/)"(-m)] 



(12) 



(13) 
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It is seen that T is the free energy of the system and that 
m is the magnetization: 



1 



(14) 



Due to the symmetry of the model, Ea. (|12|l admits two 
equivalent solutions with ±to. The choice between them 
is done spontaneously. 

The energy U and the entropies of the spins Ss and the 
occupations S are calculated from Eas. l3ll0|l 



U =- 



In ^(m) 



z(m) 



0"(/cm)ln0(fcm), (16) 
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z{m) ^-^k=±i 



(/)"(fcm) ln(j){km) - ^— ,(17) 



where z(m) = X)fc=±i'^"(^"^)- "^^^ terms in RHS of 
(|15|l are equal, respectively, to the spin-spin interaction 
energy 2^((Ei ^i^if) and aM , where = is the 

average number of particles. 

In the derivation of p 5m 7(1 we assumed that the sys- 
tem is in the ferromagnetic phase: m 7^ 0. For the para- 
magnetic phase m = one has from (|13|l 



e"'3 + 1 ' 



In 2, S 



al3 



e"/5 + 1 



ln(l 
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The phase diagram of the model is constructed in terms 
of three adimensional parameters: 



a = a/ J, 



Ts/J, e = T/j. 



(18) 



Let us start with few particular cases. For — > we get 
from ifT^ (assuming m > 0) 

2 

m = (1 + e("-™)/e.)-i(i _ _ (20) 

where C,{x) is approximated by the step function: C,(x) = 
1(0) for X > 0(x < 0). We see that m 7^ only for a < 1. 
However, only for a < 0.5, Eq. (|20|l predicts first-order 

phase transition to the ferromagnet at 9s = oi'^\ Above 
this temperature there are no particles in the system, 
TV = 0, since the energy cost for consuming a particle 
from the particle reservoir is too high (due to T — > 0). 

At 6s = 9^'^^ a finite fraction of particles is consumed from 
the particle reservoir making up the ferromagnetic phase, 
which appears as a metastable state at a higher critical 

(c) 

temperature 9* > 9s ■ There is no first-order transition 
for a > 0.5, though for 1 > a > 0.5 the ferromagnet can 
exist as a metastable phase. 

In a different limit a = we get from (|13|l : 

1' Tn Tn 771 

-^ — -9s lncosh(— ) - 01ncosh(-), (21) 

Tfi Tfi 

2m = tanh(-)+tanh(— ), (22) 
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FIG. 1: The phase diagram on the T-Ts plane; ^ = 0.45. P 
and F refer to the paramagnet and ferromagnet, respectively. 
Bold line: first-order phase transitions. Normal line: second- 
order phase transitions. For T — > 00 this line monotonically 
saturates at Tg = 0.25 J, as seen from 12311 . Dashed line: the 
instability line of the paramagnet. On the AC segment the 
latent heat is anomalous. 



Now there are only second-order phase transitions for the 
same reason as for the usual (Curie- Weiss) mean-field 
ferromagnet: d^F can never turn to zero for m 7^ 0. 
The phase diagram on the 9s — 9 plane amounts to the 
ferromagnet (paramagnet) located below (above) the line 

9-^ + 9-^ ^ 4. (23) 

The same Eas. (|21l I22II without the terms containing 9 
apply for 6* ^ 00. Here fi fluctuate so strongly that the 
influence of a disappears. 

A necessary condition for second-order transition is the 
local instability of the paramagnet: 9,^„J^|m=o = 0, i.e. 



= [1 



(2 



„a/e_^-l)]-l 



(24) 



We could continue this reasoning and develop the 
Ginzburg-Landau expansion, but here it is easier to study 
(|12l 1131) directly. For a > 0.5 there is only the second- 
order phase transition to ferromagnet at temperatures 
given by (|24|l . Ea. H14|) shows that m can increase by 
lining up the spins with or without increasing the occu- 
pations. For a > 0.5 the second way works, since it is 
too costly to absorb particles from the reservoir. Indeed, 
slightly below the transition the difference in the average 
number of particles is small: AJ\f cx tit? ^ m; see (|15|l . 

For 0.38 < a < 0.5 the phase diagram is of the type 
presented in Fig^ Now there are first-order phase tran- 
sitions from the paramagnet to the ferromagnet. Recall 
Eas. ljiyi I20|l for imderstanding their mechanism. The 
transition (phase coexistence) line is found by solving 
(|14|l , looking for the deepest minima of J-(m) , and requir- 
ing the continuous change of the free energy: !F{m) = 
JF(0). The ferromagnet corresponds to T{m) < !F{0) 
(recall the saddle-point method.) Other quantities of in- 
terest change by jump. The ferromagnet first appears 
as a metastable phase above the first-order transition 
line, while the paramagnet survives as metastable till the 
dashed line on Fig^ Thus, in the vicinity of the first- 
order transition line, both phases are locally stable. This 
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FIG. 2: The same as in Fig.JTJ, but now a = 0.245 J. 

bistability implies hysteresis and memory: when chang- 
ing the temperature the final state of the system (para 
or ferromagnet) depends on its initial state. Since both 
m and N increase in the ferromagnet, m jumps to a non- 
zero value not simply due to lining up the existing spins, 
but also due to absorbtion of additional particles from 
the reservoir. 

On the AC segment the system moves from the para- 
magnet to the ferromagnet upon increasing T or decreas- 
ing Tg. Thus, 5s is smaller in the ferromagnet, while S is 
larger there: SSg < and 5S > 0; see (O- This is neces- 
sary for the existence of the anomalous latent heat effect; 
see and below. On the AB segment the latent heat is 
anomalous: when decreasing Tg for a fixed T the system 
moves towards the ferromagnet and its internal energy 
increases. This is due to a positive energy brought about 
by the particles coming from the reservoir; see pSI) . The 
maximal magnitude of the effect reached in the middle of 
AC segment is 6U ^ J/20. Recall that the transition to 
the ferromagnet is induced by the tendency of the spin 
interaction energy to decrease. The latent SU is a con- 
sequence of this tendency that for 5U > overcomes its 
cause. On the whole first-order transition line the to- 
tal entropy iStot = S + Ss, given by IjlGI 117(1 . jumps by 
a negative amount: expectedly the ferromagnet is more 
ordered than the paramagnet. 

At the point B the latent heat is zero, while the BC 
segment shows an anomalous latent heat in a different 



scenario: upon decreasing T the system goes to the para- 
magnet, where the energy is higher. The main difference 
between the two scenarios is that now the low-T phase 
(paramagnet) has higher total entropy iStot^ a first-order 
transition to a higher entropy phase is induced by de- 
creasing the temperature T. This is impossible in equi- 
librium; see Eq.(|5J|. 

An example of the phase diagram for a < 0.38 
is presented in Fig|21 The qualitative shape of the 
ferromagnet-paramagnet boundary is given by H23|) . The 
transitions to the ferromagnet are always induced when 
decreasing either temperature: there is no anomalous 
latent heat here. During the first-order transitions the 
particle number still increases: SJ\f > 0, but now this 
brings to decreasing entropy SS < 0. In the present low- 
a regime the occupations are slaved by the spins so that 
the incoming particles arrange their entropy in the way 
dictated by the spin-spin interaction energy. For j ^ 0, 
the first-order transitions gradually disappear. 

In conclusion, we found first-order phase-transitions 
with an unexpected sign of the latent heat: it is posi- 
tive when going from high temperature (high entropy) 
phase to the low temperature (low entropy) one. Two 
conditions are necessary for this: global non-equilibrium 
(two different temperatures), though the local equilib- 
rium is kept and leads to a slightly generalized thermo- 
dynamics, and conflicting interactions. The effect was 
displayed on the model of lattice gas whose particles carry 
ferromagnetically interacting Ising spins. Motivating by 
NMR/ESR physics we assumed that the spins equilibrate 
at a temperature different from the lattice one. For this 
model there is a conflict between attractive interaction 
facilitating the income of particles from the reservoir and 
the chemical potential that makes this income costly. 

The effect seems to be generic, since we found it in 
many other situations, e.g., in the present model, but 
with fast spins and slow occupations, models with multi- 
spin interaction, the Ising spin-link ferromagnet, etc. 
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